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Abstract. Stars twinkle because their light goes through the atmosphere. The same phenomenon is expected when 
the light of extra-galactic stars goes through a Galactic - disk or halo - refractive medium. Because of the large 
distances involved here, the length and time scales of the optical intensity fluctuations resulting from the wave 
distortions are accessible to the current technology. In this paper, we discuss the different possible scintillation 
regimes and we focus on the so-called strong diffractive regime that is likely to produce large intensity contrasts. 
The critical relationship between the source angular size and the intensity contrast in optical wavelengths is 
also discussed in detail. We propose to monitor small extra-galactic stars every ~ 10 s to search for intensity 
scintillation produced by molecular hydrogen clouds. We discuss means to discriminate such hidden matter signal 
from the foreground effects on light propagation. Appropriate observation of the scintillation process described 
here should allow one to detect column density stochastic variations in Galactic molecular clouds of order of 



3 X 10"^ g/cm^ that is 10^^ molecules/cm^ 



per ~ 10 000 km transverse distance. 
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1. Introduction: The Galactic hidden matter 
problem and the baryons 

The study of rotation curves of spiral galaxies has led to 
the hypothesis of massive extended halos of "dark matter" 
(e.g. Primack et al. 1988). From the big-bang nucleosyn- 
thesis theory and the measured primordial abundances of 
the light elements, it is also established that baryons in 
the Universe are at least ten times more abundant than 
in the visible matter (stars, dust and gas) (Olive 2000). 
This deficit of baryonic matter in the Universe is ap- 
proximately equal to the deficit of gravitational matter 
in our Galaxy. This coincidence was one of the major 
motivations for the first microlensing searches for bary- 
onic hidden matter under the form of Massive Compact 
Halo Objects (MACHOs). Considering the results of 
these searches (Lasserre et al. 2000; Afonso et al. 2003; 
Alcock et al. 2000), it seems that the only constituent that 
could contribute quite significantly to the Galactic bary- 
onic hidden matter is the cool molecular hydrogen (H2). 
Indeed, this material is very difficult to detect, due to 
the symmetry of the H2 molecule that cancels the dipo- 
lar electric transitions, the only ones to be excited in a 
cold medium (~ 10 K). The perspectives for detecting 
cold H2 in outer Galactic disks have been reviewed by 
Combes & Pfenniger (1997). Most of the techniques cur- 



rently used to estimate H2 Galactic contribution are indi- 
rect detections using tracers such as CO molecule or dust, 
that imply specific hypotheses. A possible direct method 
using the ultra-fine transitions of the nuclear spins (from 
parallel to anti-parallel) concerns radio emission or ab- 
sorption (A — 0.5 km and 5.5 km) that is impossible to de- 
tect in the terrestrial environment. Detection of Galactic 
H2 clouds in front of quasars is not easily feasible, be- 
cause the absorption lines have wavelengths shorter than 
110 nm. As a summary, the hypothesis of a leading contri- 
bution to the halo mass due to cold molecular hydrogen 
is not yet ruled out by any of the currently used meth- 
ods. A hierarchical structure for cold molecular hydro- 
gen has been suggested by Pfenniger & Combes (1994) 
to fill the Galactic thick disk and De Paolis et al. (1995 
and 1998) have considered cold molecular hydrogen clouds 
as possible candidates for the Galactic halo dark matter. 
According to Pfenniger & Combes (1994) model, the gas 
could form "clumpuscules" of 10 AU size at the smallest 



scale, with a column density of 10 



,24-25 



cm and a sur- 
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face filling factor less than 1%. We propose to search for 
such cool molecular clouds in the thick disk and in the 
halo through their diffraction and refraction effects on the 
propagation of remote stars light. The method that will be 
discussed in this paper is of more general use and could 
be used to detect other types of transparent structures. 
Nevertheless, our immediate aim is to demonstrate its fea- 
sibility and its sensitivity to gaseous structures that are 



2 



M. Moniez: Does transparent hidden matter generate optical scintillation? 



considered as viable candidates for the Galactic hidden 
matter. 

2. Detection mode of extended H2 clouds: 
principle 

Due to index refraction effects, an inhomogeneous trans- 
parent (gaseous) medium distorts the wave-fronts of in- 
cident olcctro-magnctic waves. When observing a remote 
source located behind a gaseous structure (hereafter called 
screen), the luminous amplitude results from the propa- 
gation of the distorted wave- front. Diffraction theory pre- 
dicts that interference patterns should form and also pos- 
sibly longer scale refraction effects (prism and lens-like). 
In the following, we will study the visibility of the wave- 
front distortions through the intensity contrasts produced 
on Earth. We will particularly discuss the critical aspect 
of the spatial coherence and show that only the light of 
faint stars (remote and small) can provide a detectable 
signal of intensity variations due to a refractive structure. 

Before developing the scintillation mechanisms, we 
stress that the effect of a refraction process is cumulative, 
in the sense that the more material the light encounters, 
the larger are the wave distortions. If there is no dust, the 
sources behind a gaseous structure are always visible pro- 
vided that the wavelength is not resonant with a transition 
process of the medium constituents. 

3. Refraction by Galactic hidden gas 

The elementary process responsible for the refraction in- 
dex effect is the polarizability of molecules. After propa- 
gation along a distance the optical path difference be- 
tween the vacuum and a medium characterized by a num- 
ber density N (number of molecules per volume unit) and 
polarizability a is 



S = 2nNla. 



(1) 



It depends only on the total column number density Nl 
of molecules and on their average polarizability a, but 
not on the details of the density distribution along the 
propagation path. In the present case, the density is so 
low (^ 10^ moleculcs/cm^) that the medium cannot be 
considered as continuous at the optical wavelength scale 
(the average distance between molecules is ~ 10 /x » A). 
Then it is not possible to consider that the polarizabil- 
ity in a small region (compared with the optical wave- 
length) results from the average of polarizabilities over 
the molecule orientations as it is usually considered in 
dense mediums^. Nevertheless, we study the result of the 
diffusion at extremely large distances from the screen 
and precisely towards the initial propagation direction, 
where the diffusion is coherent: then the collective effect 



^ As the H2 molecule is not spherical, its polarization may 
not be aligned with the electric field, and the polarization av- 
eraged over the molecule orientations is different than the in- 
dividual polaxizations. 



of successive diffusions by randomly oriented molecules 
scattered along the propagating path is the same than 
the diffusion by molecules with the average polarizability 
(Hamilton 1960). 

For a gas only made of H2 molecules, a = 0.802 x 
10-24 cm3 (CRC Handbook 1998-99)2. 

If the Galactic halo is as described by 
Caldwell & Coulson (1986) (the so-called standard spher- 
ical halo) and if it is completely made of H2, then the 
average column density from the Sun to LMC (respec- 
tively SMC and M31) is 0.0253 g/cm^ (resp. 0.0307 and 
0.0206), corresponding to Nl = 7.62 x 10^-'^ molecules per 
cm^ (resp. 9.25 x 10^^ and 6.20 x 10^^) or to a column of 
2.83 m (resp. 3.44 m and 2.31 m) of H2 under normal pres- 
sure and temperature conditions. The extra optical path 
induced by this medium with respect to the vacuum is 
6 ~ 0.4 mm ~ 768A at A = 500 nm towards LMC (resp. 
932A and 625A). The same orders of magnitude are ex- 
pected if hidden matter lies in the thick disk instead of 
the halo. 

These values are average values and we have to 
take into account the structuration of this gas. In the 
Pfenniger-Combes model, the smallest H2 structures are 
10 AU wide and have a Jupiter mass. Therefore their sur- 
face filling factor is less than 1%. This means that we 
expect concentration factors of the column density of 100 
at least (^ Nl - 10^^ cm^^) for 1% of the sky fields. For 
such a structure the average transverse gradient of optical 
path differences is of order of 



800 x A x 100 {concentration) 



IX per 10 000 km, 



5AU 

for A = 500 nm. As we will see in the next section, this fits 
with the typical optical path fluctuations that can produce 
interference patterns for a H2 cloud located in the Galactic 
disks or Galactic halo. 

4. Diffractive scintillation: Fresnel diffraction 
applied to a stellar source and a Galactic screen 

4.1. General formalism, notations 

Let ^2(2^2, 2/2) be the luminous (complex) amplitude pro- 
duced in the source plane (see Fig. 1). For a monochro- 
matic (wavelength A) point-like source of amplitude A2 
located at (.^2,2/2) in the source plane, omitting the time 
periodic factor e"**^*, the amplitude on the screen before 
diffusion is given by the spherical wave equation: 

Ai(xi,2/i) = A2e^'='--/ri2 A2e^'=''"/0i, (2) 

where k is the wave vector (fc = 27r/A), and 

ri2 = y^?T(^^^"^^)2T(ir^^^ (3) 

= , /l + (^l^)2 + (^l^)2 
V Zi Z\ 



^ For a H2/He mixing with 24% He by mass - corresponding 
to the primordial abundances - < a >= 0.720x lO"^* cm"*, and 
subsequent calculations should be renormalized accordingly. 
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Fig. 1. Notations: The source is located in the (2:2, 2/2) plane, 
the screen contains the diffusive structure, and the observer is 
located in the {xo,yo) plane. Ai{xi,yi) and Ai{xi,yi) are the 
amplitudes before and after screen crossing. 



can be approximated by zi in the denominator because 
we will always be in the situation where zi S> xi and yi. 
The effect of the screen on the wave propagation can be 
represented by a phase delay that depends only on xi and 

yi- 



A'iixi,yi) = 



^ik5{xi,yi) 



Zl 



(4) 



where d{xi, yi) is the extra optical path due to the screen. 
The amplitude after subsequent propagation in vacuum is 
given by the Huygens-Fresnel diffraction principle: 



Ao{xo,yo) 
where 



+00 



^i(a;i,yi)-T — dxidyi , (5) 

l\Zo 



roi = zo\ 1 + ( V + ( 



zo 



(6) 



4.2. Approximations 



The stationary phase approximation states that the main 
contribution to integral (5) comes from the {xi,yi) do- 
main where the phase term of the integrand does not os- 
cillate too fast. In this domain, the Fresnel approximation 
is vahd. It consists in keeping only the first order develop- 
ment of the square root in the expression of rpi : 



''01 - Zq 



2 Zq 2 Zo 



(7) 



This approximation is in principle only valid if the next 
order terms are negligible. But as the integrand oscillates 
very fast as soon as {xq — xi)'^ + {yo — yi)^ > Azq, the 
contribution of the integral vanishes before the Fresnel 
approximation fails. This is why we can keep the infinite 
limits of the integration domain, without the "artificial" 
need of a pupiP . Within the framework of Fresnel approx- 
imation, the amplitude in the observer's plane is 

^o(a;o,2/o)= (8) 

^ This property has been checked in the simple example con- 
sidered below. 



gifczo 

iXzQ 
2iTrRj, 



4-oc 



A[{xi,yi)e 



(xQ-xi) +(!<o-!<i) 



dxidyi , 



where Rp — ^ z^/k — ^\zqI2ti is the Fresnel radius. 
Rf is of order of 1500 km to 15 000 km at A = 500 nm, 
for a screen located between 1 kpc to 100 kpc. This length 
scale characterizes the (xi, y\) domain that contributes to 
the integral (a few Fresnel radii). Then it is clear that 
the Fraunhofer approximation - that assumes the con- 
tribution of [x\ + yl)/2Rp to the phase term develop- 
ment as a constant over the pupil - is not valid here. In 
return, as we will consider sources located much further 
from the observer than the screen, we can approximate the 
spherical incident wave to a plane wave. For a point-like 
source located at {x2,y2) = (0,0), we therefore assume 
that Ai{xi,yi) = A is a, real constant just before the 
screen. The amplitude in the observation plane is then 



^o{xn,yo) = 
2iTrR% 



(9) 



gi/c5(ani,j/i)g 



1) +(yo-yi 



dxidyi 



The resulting amplitude is affected by strong interferences 
in the observer's plane (the so-called speckle) if 6{xi,yi) 
varies stochastically of order of A within the Fresnel do- 
main. This is precisely the same order of magnitude as the 
average gradient that characterizes the cool H2 structures. 

The above expression corresponds to a point-like 
monochromatic source. We will now consider a simple 
configuration to discuss spatial and time coherence effects 
that severely limit the visibility of diffraction fringes. 



5. DifFractive scintillation: a simple, pedagogical 
(and realistic) screen example 

Fresnel formalism cannot provide a simple procedure to 
take into account the source size. In this section, we an- 
alyze the effect of the source size in the case of a simple 
phase screen. Let us assume that the screen is a step of op- 
tical path S parallel to the j/i axis, described by a Heaviside 
distribution in the (xi,yi) plane 6{xi,yi) = S x H(xi). 
This case is realistic since, at the Fresnel scale, the edge 
of a gaseous structure can be considered as a straight 
line that divides the plane into two regions. The step ap- 
proaches the effect of a "strong" local gradient of the op- 
tical path. The case of a ramp instead of a step has also 
been examined, but we will only show here an example of 
diffraction pattern, because the calculation of the source 
size effect is more complicated in this case. A more com- 
plete discussion is proposed below, that concerns a wide 
domain of screen models with stochastic optical path vari- 
ations. Nevertheless, this example will be our guide for 
feasibility studies. 
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5.1. Point-like, monochromatic source 

The integral (9) that corresponds to a monochromatic 
point-source can easily be separated into a product of two 
integrals. The integral along y can be estimated by notic- 
ing that in absence of screen there should be no effect on 
the propagation. We thus get the following relation 



+ 00 .- (3=0-3=1) 



dxi 



= 2iTTRi 



(10) 



Splitting the integral along x into two parts (for xi < 
and xi > 0) we get the expression 



Aaixa,yQ) 



(11) 



Ae 



+ 00 



Ak5 



l]e 



(a=Q-a:i) 



dxi 



It follows that the amplitude in the observer's plane can 
be expressed by: 



Ao(xo,yo)=^e''""x 



(12) 



1 + + S{Xo) + C(Xo) - z(5(Xo) - C(Xo))] 



where is the reduced variable, defined by 

-^0 = Xa/iy'TTRF), 

S and C are the Fresnel integrals 



X 



Six). 

C{X) = ( 
Jo 



sm at, 

2 

cos 



The intensity in the observer's plane is 

Io{xo,yo) = AQ{xo,yo) x A'^{xo,yo) = x zo(^o) 

where 

zo(Xo) = l-{SiXn)-CiXo))sm{k5) + 



(13) 

(14) 
(15) 

(16) 
(17) 



s{Xor+c{Xor-- 



[l - cos(fc(5)] 



Fig. 2 displays the variations of this intensity in the ob- 
server's plane for S — A/4, and shows the contrast of the 
diffraction pattern as a function of the step size S. The 
inter-fringe is -in a natural way- close to the length scale 
defined by y/nRp. Fig. 3 also shows the diffraction pat- 
tern produced by a prism whose edge is the yi axis, with 
an optical path variation S — A/2 per one unit of ^/^TRp. 
Such a gradient is comparable with the average gradient 
expected from the width variation of 10 AU gas structures 
in the Galactic thick disk. One should notice here that 
diffraction patterns take place if the second derivative of 
the optical path is different from zero. Discontinuity (of 
the optical path or of its derivatives, as in our examples) 
is not necessary to get such patterns. 




- 2 




1 r ' _ ' ' ' 
fringe 
0.8 contrast 

0.6 

0.4 

0.2 
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0.2 0.3 
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0.4 0.5 



Fig. 2. Top: Diffraction pattern produced in the observer's 
plane, perpendicularly to a step of optical path i5 = A/4. The 
Xq — axis origin is the intercept of the source-step line with 
the observer's plane. Xo = 1 corresponds to xo = ^/ttRf. 
Middle: Diffraction pattern as a function of the step size 5 (in 
depth, from 0. to 0.5 in A units). 

Bottom: Contrast of the diffraction pattern as a function of the 
step size S. 



5.2. Disc-source. Spatial coherence 

A simple geometrical construction suggests that the 
diffraction pattern produced in the observer's plane by 
a point source located at {x2,y2) is the same as the one 
produced by a point source at the origin, but translated 
by i~X2Za/ zi, '-y2Zo/ zi). Actually, it is easy to check 
that, within the Fresnel approximation - using expres- 
sions (3), (6), (7) - , the difference between the phase 
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Fig. 3. Diffraction pattern produced perpendicularly to a 
prism of optical path with 5 = Xo x A/2 for Xo > 0. The 

Xq— axis origin is the intercept of the source-corner line with 
the observer's plane. Xq = 1 corresponds to xo = ^/wRf- 



fcroi + fcri2 + S{xi,yi) along the optical path defined by 
{x2, 2/2), (■'£'! ! JJi) and (xo, yo) and the phase along the path 
defined by (0,0), and {xq + X2Z0/ zi,yo + y2Z(-)/ zi) 

is independent of the point {xi,yi). This constant phase 
difference term can be factorized in integral (8), and it fol- 
lows that it is the only difference bctwcc^n tlic^ amplitude 
diflFracted at {xq, yo) from a source located at (x2, 2/2) and 
the one diffracted at {xq + X2Zo/zi,yo + y2Zo/zi) from a 
source at the origin. The intensity of the diffracted wave 
from a point-source located at {x2,y2) is then given by 



I()ixQ,yo)[source at {x2,y2)] 



(18) 



-X2^)[source at (0,0)]. 



An astrophysical source can (almost always) be considered 
as a uniform disk (radius Vg) of incoherent point-source 
elements. We then have to integrate the intensities from 
the source elements to get the resulting diffraction pattern. 
As elements aligned parallel to the y-axis give the same 
diffraction pattern, the total intensity is given by 



iixo,yo) = 

irri 



Xo 



-X2^ 
^ zi 



dx2{l9) 



where / is the total intensity emitted by the object on 
the screen and io is given by expression (17). Using the 
reduced variables 



Xn = 



and Xo 



■ 



£0 
Zl 



leads to 



I{xo,yo) = 

where 
Rs = 



21 



Rs 



Rs 



X2 

Rs 



zo(Xo + X2)dX2 ,(20) 



(21) 



is the reduced source radius which is close to the pro- 
jected source radius onto the screen plane when zi ^ zq, 
expressed in units of ^JttRf- It appears that the diffraction 
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Fig. 4. Diffraction pattern produced perpendicularly to a step 

of optical path of 5 = A/2, for extended disk-sources of reduced 
radius Rs — 1 (top), 2 (middle) and 4 (bottom). The size of 
the central depression corresponds to the reduced source size. 
Notice that the vertical scales are very different because the 
contrast is much lower for Rs = 4. 



patterns from the different parts of an extended source will 
be washed out if Rs 3> 1 or, equivalently, if the angular 
radius of the source 9s is much larger than the angular 
Presnel radius Op- This is why there is a possibility to 
observe a reasonably contrasted pattern only for extra- 
galactic sources, as shown in table 1. In this table, we 
have focused our attention on the hottest stellar classes 
that corresponds to the smallest stellar radii. 

As an illustration. Fig. 4 shows the diffraction patterns 
for a source of reduced radius equal to 1, 2 and 4. Fig. 5 
gives the maximiim contrast of the pattern as a function 
of Rs for two values of the step size S. This shows the very 
strict limitations on the angular size of the target object 
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Table 1. Examples of sources at different distances. 

- For LMC/SMC and M31 we have chosen two series of sources with comparable apparent luminosities (My ~ 18.5 and ~ 20.5). 

- The size of a SNIa envelope is estimated assuming lOOOOkm.s"' expansion rate during 20 days (Filippenko 1997). 

- The angular and luminosity distances of the Einstein cross multi-image quasar are estimated with the cosmological parameters 
of Bennett et al. (2003). The upper size of the continuum source is constrained by microlensing studies. The absolute magnitude 
given here is only indicative, as it does not take into account the magnification due to gravitational lensing. 







Distance 


Source or 


Absolute 


Apparent 




Angular 


size 9s 


Source 


Distance 


modulus 


stellar type 


mag. V 


mag. V 


Size 


rad 




Nearby star 


10 pc 


0. 


M5V 


12.3 


12.3 


O.27r0 


6.1.10""' 


1.3.10"* 


Galactic star 


8kpc 


14.5 


KOV 


5.9 


20.4 


O.85r0 


2.4.10"'^ 


5.0.10"^ 


LMC star 


55kpc 


18.7 


B8V 


-0.25 


18.45 


3.0ro 


1.2.10"'^ 


2.5.10"' 








A5V 


1.95 


20.65 


1.7r0 


7.10"'=' 


1.4.10"^ 


SMC star 


58kpc 


18.8 


B8V 


-0.25 


18.55 


3.Or0 


1.1.10"'^ 


2.3.10"' 








A5V 


1.95 


20.75 


1.7r0 


7.10"'-'' 


1.4.10"^ 


M31 star 


725 kpc 


24.3 


05V 


-5.7 


18.6 


12.ro 


3.7.10"''' 


7.6.10"* 








BOV 


-4.0 


20.3 


7.4r0 


2.3.10"" 


4.7.10"^ 


SNIa@max 


0.9 Gpc 


39.8 


SNIa 


-19.2 


20.6 


~ 2.10'" km 


7.2.10"'=* 


1.5.10"' 


at z=0.2 


















Einstein cross A 


1.75 Gpc 


45.5 


quasar 


28. 


~ 17. 


< 6.10'" km^ 


< 1.1.10""' 


< 2.3.10"'^ 


z=1.695 












< 3.10" km** 


< 5.6.10"'^ 


< 1.1.10"^ 



" Wyithe et al. (2000). 
Atsunori (2001). 



Table 2. Examples of screen positions, corresponding Fresnel sizes, and fringing time scales at A = 500 nm. The typical relative 
velocity of the nearby stars is taken from the dispersions published by Cox (Allen) (2000) 







solar 


solar 


Galactic 


Galactic 


Galactic 


Screen type 


atmosphere 


system 


neighbourhood 


thin disk 


thick disk 


halo 


Distance 


10 km 


lAU 


10 pc 


300 pc 


1 kpc 


10 kpc 


Fresnel size 


2.8 cm 


109 m 


157 km 


860 km 


1570 km 


5000 km 


Angular size Of (rad.) 


2.8.10"" 


7.3.10"'" 


5.1. 10"'=* 


9.3.10"'** 


5.1.10"'* 


1.6.10"'* 


Of in " 


0.6 


1.5.10"* 


1.0.10"^ 


1.9.10"* 


1.0.10"* 


3.3.10"^^ 


relative; speed 


Im/s 


10 kni/s 


20 kni/s 


30 km/s 


40 kni/s 


200 km/s 


time scale 


0.03 s 


O.Oi s 


8s 


2iis 


40 s 


2."3s 



to get a measurable contrast: Os should not be too large 
with respect to Op- 

Examination of tables 1 and 2 shows that every star 
gives a contrasted diffraction pattern through the atmo- 
sphere if the optical path changes by the order of A at the 
cm scale. This happens for a star close to horizon, and 
explains the stronger twinkling contrast observed at this 
moment*. More interestingly, we note that only remote 
stars have a chance to give a reasonably contrasted fringe 
system on Earth through a Galactic screen (at > 300 pc). 
We will focus on two configurations: A small source lo- 
cated in a Magellanic Cloud (LMC or SMC) or a bigger 
source located in M31, with a screen located at a typical 
halo object distance (~ 10 kpc). 

Numerical calculations show that the maximum con- 
trast region has the reduced size of the source (see Fig. 
4). The maximum intensity occurs when the source pe- 
riphery is almost tangent to the step when seen from the 
observer's plane. When the step is out of the projected 

* When the stars are high in the sky, their twinkling origi- 
nates in the weak scintillation regime described below. 



source disk, the inter-fringe is ~ 1.5 x Rp, independently 
of the source radius. 



5.3. Chromatic source. Temporal coherence 

Up to now we have only considered monochromatic 
sources. The standard UBVRI filter system has passbands 
that all satisfy AA/A < 0.1. With these filters, tempo- 
ral coherence is sufficiently high to enable the formation 
of contrasted interferences between different parts of the 
wave- front, as long as the optical path difference is less 
than A^/AA ~ 10 x A. As this coherence length is large 
compared to the typical optical path differences we con- 
sider here, we can ignore the temporal coherence aspects 
in the following. 

The inter- fringe scales with Rp, i.e. with V^, then 
AA/A < 0.1 => ARp/Rp < 0.05. The contribution of the 
wavelength dispersion to the fringe jamming is then usu- 
ally much smaller than the contribution due to the source 
extension in the configurations studied here. 
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01 2345678 



Reduced source radius Rg 

Fig. 5. Fringe contrast as a function of the reduced star radius 
Rs (~ ds/\/nOF when the source is much farther from the 
observer than the screen), for 2 values of the step size: 5 = \/A 
(thin line) and 5 = A/2 (thick line). For a large source radius, 
the contrast varies like 1/Rs. 



5.4. The observable phenomenon, summary 

Our study shows that an interference pattern with intcr- 
fringe of ~ Rp (100 — 10 000 km) is expected on Earth 
when the line of sight of a sufficiently small astrophysi- 
cal source (such as a remote star) crosses the edge of a 
structure that changes the optical path by a significant 
fraction of A. Such structures move with respect to the 
line of sight with typical velocities given in tabic 2. The 
observer's plane is then illuminated by an interference 
pattern moving with the same velocity (remember that 
the source is much farther than the screen). The shape 
of the interference pattern can also evolve, due to ran- 
dom turbulence in the scattering medium. We will base 
the present stiidy on the assumption that the scintillation 
is mainly due to pattern motion rather than the pattern 
instability, as it is usually the case in radioastronomy ob- 
servations (Lyne & Graham-Smith 1998). The contrast of 
the pattern is critically limited by the angular size of the 
source; in the configurations proposed here, we expect a 
typical contrast ranging between 1% and 10%. The time 
scale tscint of the intensity fluctuations is tsdnt = RfIVti 
where Vt is the transverse velocity of the structure. It is 
of order of 10 — 40 s. As discussed just before, the inter- 
fringe scales with -s/A; therefore, one expects a significant 
difference in the inter-fringe and in the time scale tsdnt 
between the red side of the optical spectrum and the blue 
side. This property might be used to sign the diffraction 



phenomenon at the Rp natural scale. We will see below 
that the A dependence is different if the screen produces 
optical path stochastic fluctuations at a length scale much 
smaller than Rp. 



6. General discussion of the scintillation regimes 

Up to now, we have considered variations of the screen op- 
tical thickness at a scale close to Rp. In the screen plane, 
Rp can be considered as a coherence domain, in the sense 
that - as mentioned earlier - integral (9) is dominated by 
the contribution of a region characterized by the Fresnel 
radius. This means that only the patch of the wavefront 
within a few Presnel radii contributes coherently to the 
integral. The details of the optical path variations inside 
this domain then drive the diffraction pattern, leading to 
the (small scale) diffractive scintillation. If there is a large 
scale structure (^ Rp) in the optical path variations, a 
succession of focusing and defocusing configurations oc- 
curs, where the first phase term k5{xi,y{) in integral (9) 
can partially compensate or enhance the second phase 
term; the coherence domain that contributes to the in- 
tegral can then become larger or smaller than Rp.ln the 
case of large optical path variations, several coherence do- 
mains can also converge and the average luminosity can 
be estimated through the geometrical optics approxima- 
tion. The intensity then depends on the local focal length 
produced by the optical path variations on the screen. 
Intensity variations - called refractive scintillation - arise 
at a length scale larger than Rp. 

Intensity scintillation is well known in observations 
of compact radio sources. Indeed scintillation of radio- 
pulsars is used to study the nearby interstellar and so- 
lar system media. The large wavelength and the very 
small size of the sources allow one a good control of 
this technique. In this section, we will summarize the 
conditions and characteristics for the different scintilla- 
tion regimes that occur in pulsar radio-observations, and 
adapt them to star optical observations. Extensive liter- 
ature exists on pulsar physics. Reviews can be found in 
Lync & Graham-Smith (1998) and Narayan (1992); the 
scaling laws described below are taken and adapted for 
the latest reference. 

The distinction between the regimes is based on the 
relative values of the length scale of the optical path fluc- 
tuations and of the Fresnel scale Rp. The screen is char- 
acterized by its diffractive length scale Rdiff, defined as 
the separation in the screen plane for which the root mean 
square of the optical path difference 6{xi ,y[) — 6 {xi ,yi) is 
A/27r (see Fig. 6). As the patch of the wavefront contribut- 
ing to integral (9) has a size of order of Rp in the screen 
plane, we can consider two very different situations since 
integral (9) has a completely different behaviour whether 
Rdiff > Rf or Rdiff < Rf- 
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Strong diffusion 



Fig. 6. The two scintillation regimes: sections of the wavefront 
after screen crossing. 

Top panel: Rdiff 3> -Rf- The weakly distorted wavefront pro- 
duces a weak scintillation at scale Rf in the observer's plane. 
Bottom panel: Rdiff ^ Rf- The strongly distorted wavefront 
produces strong scintillation at scales Rdi / / (diffractive mode) 
and Rref (refractive mode) in the observer's plane. 



6.1. Weak scintillation regime 

If Rdiff ^ Rf, then the random fluctuations of the op- 
tical path are small compared to the wavelength within 
the Fresnel zone. The screen is then qualified as weakly 
diffusive. Only weak perturbations of the wavefront are 
expected, essentially from the focusing/defocusing effects 
of the optical path fluctuations within this Fresnel zone. 
In this regime, only the wavefront patch enclosed in the 
Fresnel radius plays a significant role. Indeed as variations 
of 6{xi,yi) are small over the Fresnel domain, d{xi,yi) 
cannot significantly change the total phase in integral (9) 
when (xi — xa,yi — yo) is close to the Fresnel circle. 

An alternative point of view is to consider that the 
emerging wave is a combination of the initial plane wave 
plus weak amplitude perturbing waves with a random dis- 
tribution of wave normals. For the observer, a point source 
located behind the screen would appear like a main spot, 
surrounded by low luminosity spots. The resulting inten- 
sity variations when the screen moves from one Fresnel- 
size zone to the next are weak, and their typical time scale 
is Rf/Vt, where Vr is the transverse velocity of the screen 



with respect to the observer-source line of sight. Table 3 
gives the conditions to observe this scintillation regime 
and its characteristics. 

6.2. Strong scintillation regimes 

If Rdiff Rf, then the random fluctuations of the opti- 
cal path produce strong random phase changes within the 
Fresnel zone. The screen is considered as a strongly dif- 
fusive filter because the wavefront is strongly perturbed. 
From our alternative point of view, one can consider that 
the contribution of the initial plane is almost completely 
redistributed after crossing the screen. The point source 
seen by the observer through the screen would now ap- 
pear like scattered spots with random intensities. In this 
situation, the intensity is expected to strongly vary at 
the sub-Fresnel scale. Moreover, large length-scale struc- 
tures in the optical path variations produce average fo- 
cusing/defocusing effects over domains much larger than 
the Fresnel size. Two scintillation modes are then to be 
distinguished: 

— Diffractive scintillation 

As the optical path variations are larger than A within 
the Fresnel radius, one expects a strongly contrasted 
fringe system (for a point-like source). By contrast 
with the simple case studied in Section 5, where no 
explicit length scale other than Rf was introduced, 
there is now a different scale of optical path varia- 
tions 6(xi,yi) characterized by Rdiff - Then, when the 
screen is shifted from a position to another located at 
distance Rdiff, the values of integral (9) are not corre- 
lated between the two positions. The resulting inten- 
sity variations when the screen moves are large, and 
their time scale is Rdiff/Vr- 

As discussed in Section 5.2, diffractive scintillation is 
critically dependent on the source size. In contrast, the 
time scale does not depend on the wavelength A, but 
the positions of the maxima/minima do, thus limiting 
the temporal coherence of the fringes. When A changes 
by AA, the fringes produced by Rdiff -size structures 
are displaced by ^ zqAX/ Rdif f - The fringe systems 
are decorrelated when this displacement is ^ Rdiff, 
i.e. when AA/A - {Rdiff /RfY (see Table 3). 

— Refractive scintillation 

This effect results from the strong stochastic variations 
of the optical path produced by large scale structures 
of the screen. Let Rref be this additional length scale. 
In contrast to the case of weak scintillation, the large 
variations of the optical path can considerably change 
the values of the integrand phase of integral (9) at 
the limits of the Fresnel domain. Focusing/defocusing 
effects are not restricted to the Fresnel zone any more, 
and light from "other" Fresnel domains can be focused 
if the optical path variations are large enough. Such 
effects now involve regions of typical size defined by 
Rref, and their time scale is Rref/Vr. 
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Fig. 7. Apparent stellar angular radius 6s for a solar-like star 
{va = tq) and for a star with Va = lOr© as a function of 
the distance. The lower line gives the Fresnel radius _Rf as a 
function of the screen's distance, illustrating the fact that the 
source has to be much farther than the screen to give contrasted 
diffractive scintillation. 



Source size is much less critical for the visibility of 
this mode. The time scale also does not depend on the 
wavelength and the illumination patterns are coherent 
within a wide passband (see Table 3). 

The strong scattering diffraction mode is the most 
promising in terms of contrast, as demonstrated in Section 
5. Fig. 7 allows one to compare the apparent stellar an- 
gular size 0s with the natural angular Fresnel scale dp 
at A = 500 nm. Fig. 8 shows the iso-contrast lines in the 
fringing system expected in the case of strong diffractive 
scintillation. Expressed in practical units, the modulation 
index as defined in table 3 is in this case 



0.071 



(22) 



500 nm J 



Rf 



zq+zi 
1 kpc 



1 pc 



The example discussed in Section 5 corresponds to a 
diffractive scintillation regime that represents a transition 
between the weak and the strong regimes. The natural 
scale Rp is here the only relevant parameter depending 
only on the screen position and wavelength. For this tran- 
sition regime, the modulation index estimated from the 
above formula is stronger than calculated in Section 5, 
because the approximations of Table 3 are not valid. For 
a BOV star in M31 {9s = 4.7.10"** ") and a screen at 
Ikpc {dp — 1.0.10"^ "), the expected contrast for diffrac- 
tive scintillation - setting Rdiff — Rp in the expression 



Fig. 8. Lines of equal diffractive scintillation modulation index 
at A = 500 nm in the source/screen distances plane. 
The modulation index is 100% (thick line), 50%, 10%, 5% 
and 1% (from top to bottom lines), assuming that the prod- 



uct 



'Rd.ff' 




7's 









is unity. This happens only for very small 

sources, and when Rdiff is not too small compared to -Rf- 
Otherwise, the modulation index has to be multiplied by this 
product. Typical distances of disks and halo clouds are indi- 
cated, as well as the distances of the two most promising tar- 
gets. 



of Table 3 - is m.^mt ~ {O^ff/Os) - (^f/^s) = 0.21, to 
be compared with the contrast of ~ 0.1 expected from a 
step of A/2 in the optical path (Fig. 5). 



6.3. Extreme scintillation 

For completeness, we should also mention an extreme 
scintillation regime, due to caustic effects, occurring 
when the observer approaches the focus of a refrac- 
tive structure. Structures like the ones predicted by 
Combes & Pfenniger (1997), with concentrations of diam- 
eter D = 10 AU that produce a total variation in the 
optical path (from edge to center) A = 50 000 x A at 
A — 500 nm, should act as optical lenses. A naive focal 
length estimate gives / ~ D'^/SA — 360 Mpc. The con- 
centration power of such a lens is clearly much too weak 
to give a measurable signal on Earth. Moreover, the time 
scale of the intensity variations should be tgcint — D/Vp — 
(10AU)/(200km.s-i) = 87 days. Other authors have dis- 
cussed in detail the possible caustic effects with poly- 
tropic models of self-gravitating gas clouds, and the con- 
sequences on background star light-curves (Draine 1998, 
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Table 3. Conditions and characteristics of the different scintillation regimes for point-like and extended sources of apparent 
angular size 6s- Rsdnt is the characteristic length scale of the illumination pattern, and tadnt is the characteristic time scale 
of the intensity variations for a screen moving at transverse speed Vr with respect to the line of sight. The modulation index 
mscint is the root mean square amplitude of the intensity scintillation, msdnt depends only on the 6di ff to Os ratio for the 
diffractivo mode. For the other scintillation modes, msdnt is given for stochastic fluctuations of the optical path length due to 
Kolmogorov turbulence in the scattering medium. 





^diff ^ Rf 
Weak scintillation 


Rdif 
Strong s 

diffractivc mode 


f < Rf 
cintillation 

refractive mode 


characteristic length 
characteristic angle 
coherence passband A A/A 


Hp 

Of 
~ 1 


iidiff 
Odiff 

{Rd^ff/Rpf 


^ref ^ i^diff 
Oref 
~ 1 


Point-source condition 


es < Of 


ds < Sdiff 


Os < Oref 


Ftscint 
tscint 


Rf 
Rf/Vt 

{RF/Rdrff)i « 1 


Rdiff 

Rdiff /Vt 
~ 1 


Rref 
Rref/Vr 
{Rd^ff/RF)i « 1 


Extended source 


es > Of 


Os > Od^ff 


Os > Oref 


Ftscint 

tscint 

'iTT'Scint 


RFies/Op) 

(RF/VTWs/eF) 

{RF/Rdiff)HOF/es)i 


Rdiff{ds/9diff) 
{Rdif f /Vt){Os /Odif f) 

Odiff/es < 1 


Rref {Os /Oref) 
(Rref /Vt) {9 s /Oref) 
{Rdif f/RF)^ {Oref /es)i 



Rafikov & Draine 2001). We will not consider this regime 
in the following. 

6.4. Regime discrimination 

The appropriate tool to connect the screen structure with 
the scintillation features will be the temporal power spec- 
trum of the intensity fluctuations. If strong diffractivc and 
refractive regimes take place together, two peaks are ex- 
pected in the power spectra. Their relative power will 
also be connected with the source size. The so-called in- 
verse problem, i.e. the determination of a density fluc- 
tuation model of the interstellar medium from scintilla- 
tion observations has been extensively discussed by radio- 
astronomers (Narayan 1988), but is beyond the scope of 
this exploratory paper. 

We will mainly discuss here the perspectives of the 
diffractivc scintillation regime. Nevertheless, we will also 
sometimes mention the refractive regime because impor- 
tant density fluctuations may occur at a much smaller 
scale than the global structure, due to gaseous turbulence 
phenomena. 

7. Possible configurations for detectable 
diffractive scintillation 

One of the main outcomes of the previous discussions is 
that only small angular sources may give a reasonably 
contrasted diffractive scintillation. We now focus on this 
diffractive process, because its modulation index can be 
large and is easy to predict with basically no hypothesis on 
the detailed screen structure. We will only assume that a 
regime characterized by Rdiff < Rf can be established, or 
that at least a transitory regime as described in Section 5 
- characterized hy Rp - can occur, for example if an inho- 
mogeneity due to a turbulent mechanism crosses the line 



of sight. Table 4 that combines data from Tables 1 and 

2 lists some configurations that should produce diffrac- 
tive scintillation and gives their characteristics. Depending 
on the diffusion strength, other regimes may take place; 
the reader may adapt the forthcoming discussion to these 
regimes, after rescaling the characteristic length, time and 
contrast according to Table 3. Optical depth may also be 
considerably larger in those other regimes because they 
are strongly dependent on the screen structure. We will 
essentially consider minimum optical depths in what fol- 
lows. 

To end the discussion about the screen configurations, 
we have to mention here that a diffusive screen located 
near the source, or in the galaxy of the source, do not 
produce scintillation, due to the size of the source. In 
such a configuration {z\ <^ zq), the approximation of 
an incident plane wave on the screen is not valid any 
more, and calculation should be redone starting with 
eq. (8). The Presnel zone is then defined by Rp = 
s/X/2n{l/zo + l/zi)-'^ ^Azi/27r instead of ^Xzq/I-k 
(see e.g. Born & Woff (1975) or Sommerfeld (1954)). The 
zq/z\ multiplier that enters the reduced radius expression 
(21) is very large; consequently, spatial coherence is com- 
pletely lost in this situation. Moreover, a diffusive screen 
dose to the source cannot increase significantly its appar- 
ent size for a simple geometrical reason: the image of a 
source is not significantly distorted by diffusion occurring 
near the emission point. We will ignore such effects in the 
following. 

8. Optical depths and event rates for diffractive 
and refractive scintillations 

In this section, we want to quantify the probability to 
observe scintillation produced by disk or halo molecular 
clouds. The 1% surface filling factor predicted in the model 
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Table 4. Configurations leading to strong diflFractive scintillation. Numbers are given for A : 
not taken into account for the magnitude estimates. 



500 nm. Interstellar absorption is 



Distance 



Rf to multiply by 



SCREEN 



atmos- 
phere 



10 km 



2.8 cm 



solar 
system 



lAU 



109 m 



solar neigh- 
bourhood 



10 pc 



157 km 



Galactic 
thin disk 



300 pc 



860 km 



Galactic 
thick disk 



Ikpc 



1570 km 



tscint to multiply by 



r A ] 


2 


Rdiff 


L 500 nm J 







0.03 s 



0.01s 



8s 



29 s 



40 s 



Optical depth Ts, 



total > 10" ' 



mscint in % to multiply by 

L 600 nn 



Rdijf 




IS. 


Rf 




r0 



100% 



32 [t4] 



2-2 [t4] 



4 1 f^^l 

L lOkpcJ 



2.2 



[lOkpc] 



SOURCE (Cox (Allen) 2000) 



Location 



Type 



V 



DIFFRACTIVE MODULATION INDEX m.cint 
(to multiply by ^A/500 nm X Rdiff/Rp) 



nearby 

star 
d = lOpc 



B8V 

A5V 
KOV 
M5V 



1.7r0 
O.85r0 
O.27r0 



-.25 
1.95 
5.9 
12.3 



Polaris 
d= 130 pc 



F7I 



lOOro 



1.97 



Galactic 
star 

d = 8 kpc 



KOIII 
G5III 
B5V 

A5V 
KOV 
M5V 



15.rQ 
lO.ro 
3.9r0 
1.7rQ 
O.85r0 
O.27r0 



15.2 
15.4 
13.3 
16.5 
20.4 
26.8 



LMC 

d = 55 kpc 



K5III 
05V 

B8V 
A5V 
KOV 



25.r0 
12.r0 

3.r0 
1.7r0 
O.85r0 



18.5 
13. 
18.5 
20.7 
24.6 



M31 
d = 725 kpc 



B5I 
05V 

BOV 
B8V 



5O.r0 
12.rQ 

7.4r0 
3.r0 



18.1 
18.6 
20.3 
24. 



<1% 

for 

every 
source 
in a 

telescope 
> Im 

(see text) 



z=0.2 
= 0.9 Gpc 



SNIa 



2.10'' 
km 



20.6 



Einstein 

cross 
= 1.75 Gpc 



Quasar 

Q2237 
-^0305 



< 6.10' 

km 

< 3.10^ 



17 



17 



10% 

20% 
40% 
- 100% 



4% 



100% 
100% 
100% 
100% 
100% 
100% 



100% 
100% 

100% 
100% 
100% 



100% 
100% 

100% 
100% 



100% 



100% 
100% 



1% 



1% 

2% 
5% 

10% 
20% 
65% 



5% 
10% 

40% 
70% 
100% 



32% 
100% 
100% 
100% 



70% 



> 45% 

>9% 



1% 
2% 
4% 
12% 



1% 
2% 
8% 
13% 
27% 



6% 
25% 
40% 
- 100% 



13% 



>8% 
>2% 



1% 

4% 
7% 
14% 



3% 
13% 

22% 
53% 



7% 



>4% 
> 1% 



of Pfenniger & Combes (1994) for gaseous structures is 
also the maximum optical depth for all the possible refrac- 
tive (weak or strong) and diffractive scintillation regimes. 
Nevertheless, we want to consider here the pessimistic case 
where the optical depth for the strong regimes is much 
smaller: in the thick disk, the typical transverse speed of 
a structure is 40km.s~"'^. Then the typical crossing time of 
such a structure should be ~ 400 days. Under the hypoth- 
esis that strong diffractive regime is expected only when 
the structure enters or leaves the line of sight^, the dura- 
tion for this regime is of order of ~ 5 minutes (time to 
cross a few fringes). Then the optical depth Tsdnt for such 



^ The validity of this hypothesis as well as the optical depth 
estimate for each scintillation regime depend on the internal 
structure model. Turbulence or any process creating filaments, 
cells, bubbles or fluffy structures should be taken into account. 



regime is at least of order of 10"'' and the average expo- 
sure needed to observe one event of ~ 5 minute duration 
is 10^ star x hr. 

If the gaseous structures belong to the Galactic halo 
instead of the thick disk, the same order of magnitude is 
also expected for the optical depth. 

We will neglect the multi-diffusion eventuality, 
given the small sky fraction occupied by the struc- 
tures. For our study, the case of a thick (longitu- 
dinally extended) screen needs a formalism adapta- 
tion only if the extension along the line of sight 
is such that Rp significantly changes within the 
structure. Several authors (Lyne & Graham-Smith 1998, 
Tatarskii & Zavorotnyi 1980) have studied very thick 
screens and found a similar behaviour to the one that pre- 
vails for thin screens. 
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9. Feasibility of observations 

The diffractive regime with Galactic hidden H2 gas (see 
the three last columns of Tables 2 and 4) can happen 
with contrasts better than 1% only if the source has an 
angular size smaller than a few ~ 10~^^rad. For a given 
surface temperature this constraint is equivalent to a con- 
straint put on the source (star, SN, quasar...) magnitude. 
Choosing objects with the highest surface luminosity in 
Table 4 indicates that the minimal magnitude of stars 
whose light is likely to undergo a few percent modulation 
index is about My = 20.5. Therefore, diffractive scintil- 
lation search needs the capability to sample every ~ 10 s 
(or faster) the luminosity of stars with My > 20.5, with 
a point-to-point precision better than a few percent. This 
performance can be achieved using a 2 meter telescope 
with a high quantum efficiency detector. 

Indeed, the optimal relative photometric precision 
A is related to the star magnitude M, the expo- 
sure time Tejjp, the seeing, the sky background magni- 
tude /z per arcsec^ and the telescope collecting surface 
Stel = ""-Ctef/^ through the following formula, taken from 
(Moniez & Perdereau 2001): 



IqM/2.\ 



^ exp — 



SteMM = 0) A2 



(23) 



0.27 -M.5W 0.25 -Fl. 13 



seeing 



1 



10- 



where $(M = 0) is the photo-electron flux collected at ob- 
servatory level by the detector in the considered passband 
per m^ collecting surface, for a M = {] star. This relation 
is valid as long as the readout noise remains negligible. 
For a telescope with two reflectors of 0.87 transmission 
each, equipped with a high quantum efficiency detector 
(~ 90%), $(Mv = 0) ~ 4.2x 10V"/s/m2 in the V-band. 
For good observation conditions {seeing = 1", n > 21.7, 
i.e. moon- less nights), the squared term of expression (23) 
is ~ 2 for stars with My = 20.5, and negligibly smaller 
for somewhat brighter stars. Then the relation simplifies 
into 



Mv—20.5 1 
A ~ 30% X 10^ — X 



1 



Vt:^ Dtei (m) 



(24) 



where Dtei is the diameter of the telescope. The photo- 
metric precision A would be 5% with 10 s exposures on 
a My = 20.5 star, with a 2 meter telescope. Neglecting 
the interstellar absorption, this is sufficient to detect the 
most contrasted (and the longest) diffractive scintillation 
events, as can be seen in Table 4. 

We wish to stress that the continuity of the monitoring 
is not a critical issue. As events last only for a few min- 
utes in our scenario, the detection efficiency will not be 
affected by data taking interruptions, contrary to the case 
of microlensing searches. The essential parameters for the 
sensitivity of a detection setup will be the integrated expo- 



sure in star x hr®, the sampling rate and the photometric 
precision. 

The variation of the diffractive scintillation pattern 
with the wavelength A should also be used as a signature of 
the process. As emphasised in Section 6.2, if Rdiff ^ Rf 
the time scale docs not depend on A, but the phase (the 
timing of extrema) changes with A. In the particular case 
where Rdiff ~ Rf, the peak of the temporal power spec- 
trum is expected to scale with \/A- For these reasons, 
multi-wavelength detection capability is highly desirable. 
The most powerful approach would be to get star spectra 
at a high time sampling rate. The equivalent in radio- 
astronomy is called the dynamic spectrum, showing the 
intensity fluctuations in a 2D (time versus wavelength) 
diagram (sec for example Gupta et al. (1994)). 

According to Table 4, M31 hot main sequence stars are 
the most promising targets, but LMC and SMC small stars 
are probably easier to distinguish from their surroundings. 

10. Foreground effects, background to the signal 

One of the challenges to extract an interstellar intensity 
scintillation signal is the disentangling of the scintillation 
due to the foreground media (atmosphere, solar system 
and nearby medium). 

10.1. Atmospheric intensity scintillation and absorption 

Atmospheric intensity scintillation should be easier to 

handle than naively expected. Dravins et al. (1997-98) re- 
port extensive studies of this phenomenon, that is con- 
nected with seeing studies. When short length scale in- 
terference patterns (a few cm) produced in the diffractive 
mode enter a large aperture telescope, the collected light 
is averaged, since the aperture acts like a low pass fllter. 
The modulation index is then much smaller than 1% for 
a telescope diameter larger than one meter. Moreover, the 
time scale of this diffractive regime is two orders of mag- 
nitude faster than the searched Galactic signal. As far as 
the weak intensity scintillation is concerned, it also in- 
duces very small modulation index niscint at 1 — 10 s time 
scale; the following formula gives the low-frequency com- 
ponent of intensity scintillation at optical wavelengths, at 
altitude site h (Dravins et al. (1997-98) part III): 



'^scintiiexp^ ^ 0.09 



D 3 



cos(Z) 



-1.75 



g 8000 m _ 



(25) 



where cos(Z)~^ is the airmass and texp is the integration 
time. For a 2 meter telescope, an airmass of 2 and an 
observatory altitude of 2400 m, we find that niscinti^ s) = 
3.10"^ and mscint(lOs) = 10~^. Thus intensity variations 

due to atmosphere will not seriously affect a diffractive 
interstellar scintillation signal of 1% amplitude or more. 

® If dead-times larger than a few minutes occur, they should 
be subtracted from the exposure because events can be con- 
tained and subsequently lost in such time intervals. This is also 
a noticeable difference from the microlensing searches. 
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Any other long time scale atmospheric effects such as 
absorption changes at the sub-minute scale (due to fast 
cirruscs for example) should be easy to remove as long as 
nearby stars are monitored together. At this time scale, 
the light of stars within a small angular distance undergo 
the same atmospheric effects at the same time, or with 
a short delay. A careful subtraction of such collective ef- 
fects should in principle considerably reduce this source of 
background. 

10.2. Intensity scintillation due to the solar systenn 
medium, and to solar neighbourhood structures 

— Interplanetary medium: from to 5.10~^pc 

As can be seen Table 4, diffractive scintillation pro- 
duced by interplanetary medium have a characteristic 
time of 10~^ s, therefore it does not disturb an hypo- 
thetic Galactic gas signal. 

The eventuality of a slow refractive regime has to be 
examined in more details: Studies of the atmospheric 
properties by (Dravins ct al. (1997-98) part I) were 
mainly done with Polaris star. This star is far enough 
(130 pc) to scintillate through interplanetary screens. 
Unfortimately, this star is big (^ lOOro) and its light 
cannot be subjected to modulations of more than 4% 
in the diffractive regime'^. Nevertheless, in the hypoth- 
esis of a refractive regime that we are discussing here, 
we know that this not very critical, and Polaris mea- 
surements should then really be a reliable indicator of 
this latter scintillation type. As the measured modula- 
tion index of this star is much smaller than 1% at any 
frequency, this gives us a good preliminary indication 
that both diffractive and refractive scintillations from 
interplanetary medium are probably negligible toward 
Polaris direction. 

— Solar neighbourhood: from 5.10"^ to 10 pc 
This conclusion also applies to solar neighbourhood 
screens (up to ~ 10 pc) for the refractive regime, but 
the size of Polaris prevents us from drawing any con- 
clusion relative to the diffractive regime. 

The solar system is embedded in a local interstellar 
cloud extending not farther than 10 — 20 pc. The av- 
erage column density of atomic hydrogen through this 
structure is less than 10^® cm~^ towards the Galactic 
center (Lallemcnt ct al. 1995, Fcrlct 1999), which is 
negligible compared with the column density of the 
molecular clouds we are looking for. If molecular or 
atomic overdensities occur in this local region typi- 
cally at 10 pc -, then according to table 4, many types 
of stars located at ~ 8 kpc, including red giants, should 
undergo a contrasted diffractive scintillation; as it can 
be deduced from this table 4, the distinctive feature of 
scintillation through more distant screens (> 300 pc) is 
that only the smallest stars are expected to scintillate. 
It follows that a strategy consisting in the simultane- 

Polaris is also a spectroscopic binary system that may 
complicate the interference system. 
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ous monitoring of many different types of stars located 
at different distances should allow one to discriminate 
effects due to solar neighbourhood gas and due to more 
distant gaseous structures. 
— The local interstellar medium: from 10 to 200 pc 
The local interstellar cloud itself is located inside a 
bubble, which is a 50 — 200pc cavity mainly containing 
hot ionized gas. Ionized hydrogen has a negligible effect 
on the visible light propagation. The remaining atomic 
plus molecular hydrogen components contribute for 
a total column density of less than lO^'^cm"^ up 
to 200 pc (Lehner et al. 2003, Lallemcnt et al. 2003). 
This is again much smaller than the column density 
due to molecular clouds. If dense atomic or molecu- 
lar structures are located at such distances or beyond, 
they start to enter our domain of interest. 

10.3. Possible sources of fake signal? 

There is no known physical process that can produce 
> 1% intrinsic intensity variations of an ordinary star at 
the minute time scale. Asterosismology involves acoustic 
modes that produce a few tens of ppm intensity varia- 
tions in the "high" frequency domain of a few minutes 
(Christensen-Dalsgaard 1998). Planet transit could also 
give a fast luminosity change, but are also expected to ex- 
hibit intensity variations smaller than 0.1%. Granularity 
of the star surface, spots or eruptions would induce much 
lower frequency intensity variations than the diffractive 
scintillation. A few categories of recurrent variable stars 
exhibit important emission variations at the minute time 
scale (Sterken & Jaschek 1996). Among them are the rare 
types UV Ceti and flaring Wolf-Rayet stars. Both types 
are easy to identify from their spectral characteristics. UV- 
Ceti are also very faint stars (absolute magnitude > 15) 
and only the closest ones could contaminate a monitoring 
sample. 

10.4. Technical limitations 

Blending of 

stars, which is common in the LMC/SMC/M31 crowded 
fields, may attenuate the measured contrast of a fringed 
pattern. Careful simulations are needed to measure the 
impact of these limitations. Adaptative optics or space 
measurements are possible ways to reduce this impact. 
Another technical difficulty is the risk of getting compli- 
cated (and fluctuating) point spread functions as a re- 
sult of the very short exposures. Such circumstances could 
make necessary the use of performant subtraction algo- 
rithm, given the very crowded environment. 

11. Observation specifications and strategy 

A careful feasibility study should be made before start- 
ing an ambitious observation program: The stability of 
light curves of nearby stars should be checked with a low 
cost telescope at the sub-minute time scale, to test the 
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control of collective intensity fluctuations due to meteo- 
rological phenomena and to explore possible limitations 
due to interplanetary gas. After such preliminary stud- 
ies, an observational program could start with the spec- 
ifications described below. These specifications are opti- 
mized for the search for diffractive scintillation in a large 
sample of stars. We showed in table 4 that supernovae 
and "small" quasars^ may also be subjected to diffrac- 
tive scintillation; but the possible smallncss of the optical 
depth (~ 10~^) and the scarcity of such sources proba- 
bly don't make them suitable targets for the first searches 
of diffractive scintillation. Moreover, one should not for- 
get that diffractive scintillation of such sources is not ex- 
pected from clouds located far beyond the Galaxy (see 
end of Sect. 7). Conversely, refractive regimes, which may 
have much larger optical depth, could be searched on any 
target. 

11.1. From Earth, with a single telescope 

— Optics 

As already mentioned, at least a 2 meter class tele- 
scope is required to search for scintillation in < 10 s 
exposures on 20 < My < 21 LMC/SMC or M31 stars. 
First order adaptative optics will lead to a better pho- 
tometric precision, but second-order adaptative op- 
tics (that corrects focusing/defocusing effects of atmo- 
sphere) is probably not really necessary since the at- 
mospheric intensity scintillation should be negligible. 

— Filters 

Monitoring in infra-red wavelength is not clearly a de- 
cisive advantage; the Presnel radius is larger, but red 
stars arc also usually larger. In contrast, monitoring 
stars with at least two different passbands should help 
to disentangle diffractive scintillation from the refrac- 
tive one. In the case of diffractive scintillation, the 
phase differences of the fringe systems between two 
different passbands will provide constraints on the ra- 
tio of the screen distance to the diffusion scale Rdiff 
(see Section 6.2). 

— Detector 

A fast, low-background readout detector is essential to 
make possible the requested high sampling rate and to 
maintain a decent useful to dead time ratio. A possi- 
ble option could to put an array of small (high quan- 
tum efficiency) frame-transfer CCDs at the focus of the 
telescope, that could be rapidly read (because they are 
small) in parallel during the exposures. 
An alternative could be to use a series of very narrow 
CCDs (short columns), with continuous slow readout 
of the lines, to collect continuous light-curves without 



* The Einstein's cross Q2237-f0305 considered in the tables 
is the quasar tiiat has the most precise size determination. It 
is also located at an optimal redshift z=1.695, corresponding 
to the maximum possible angular distance according to the 
currently published cosmological parameters, and then to a 
minimal angular size 6s- 



any dead-time. Such continuous readout option is not 
viable with a normal (not narrow) CCD, because of 
the light-curves mixing due to field crowding. 

— Targets and fields 

Imaging 1 deg^ field on isophot My = 23 mag/arcsec^ 
towards LMC/SMC/M31 allows one to monitor of or- 
der of 10^ (20 < My < 21) stars^ (Elson et al. 1997; 
Hardy 1978). It is essential to choose fields containing 
stars with radii spanning a wide interval, and also con- 
taining a sample of nearby stars. This is necessary to 
extract information about the screens through the re- 
lationship between the temporal power spectrum and 
the source size and distance. In particular, nearby stars 
should provide a control sample of non-scintillating 
light-curves. 

— Telescope time 

Exposure of a few 10^ star x hr (on small stars with 
20 < My < 21) during dark nights would provide a 
significant sensitivity to the existence of Galactic II2 
structures. This needs typically one season of obser- 
vations using a wide field telescope (1 deg^) or a few 
seasons with a standard field telescope. Nevertheless, 
as already mentioned, there is no need for consecutive 
telescope time, making such program relatively flex- 
ible. Every sequence of consecutive images - longer 
than a few minutes - will contain its own events and 
will then be autonomous. Distinct but possibly simul- 
taneous campaigns with different air-masses could cer- 
tainly improve the knowledge of atmospheric effects. 

— Complementary observations 

In case of positive detection, complementary observa- 
tions should be planned for scintillating candidate ob- 
jects. The first requirement is the identification of the 
stellar type through spectroscopy, in order to check 
that the candidate do not belong to the very specific 
cataclysmic objects mentioned in section 10.3, and to 
get an estimate of the distance and radius of the can- 
didate. Such estimates are necessary to constrain the 
maximum distance of the screen. 

In the case of a star appearing to scintillate for a long 
time, time-resolved spectroscopy or multi-wavelength 
data taking as well as multiple detection (see below) 
will allow one to check the hypothesis of diffractive or 
refractive scintillation. 

11.2. From Earth: array of telescopes 

A 2D array of telescopes, a few hundred and/or thousand 
kilometers apart, would certainly be the most powerful 
system to measure the characteristics of a diffusive screen: 
One single telescope will only provide a degenerated in- 
formation (the time scale tsdnt = Rdif f /Vt)- Sampling 
a diffraction pattern with a 2D array would allow one to 
separately measure the geometrical scale Rdiff and the 
speed of the pattern. Atmospheric effects will be decorre- 

® In more luminous fields like in the LMC bar, these stars 
would be blended by more luminous ones. 
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lated between the telescopes, as well as the interplanetary 
gas effects. Nearby gaseous structures (~ 10 pc) are ex- 
pected to produce ^150 km fringes on Earth and Galactic 
structures are expected to produce > 1000 km fringes, that 
should be easy to distinguish with a handful of synchro- 
nized telescopes sampling a few thousand kilometers wide 
pattern. The speed and direction of the diffraction pat- 
tern (drift due to the relative velocity of the screen with 
respect to the line of sight) . as well as the pattern's vari- 
ations (due to the dynamics of the scattering medium), 
should also be measurable through the analysis of time 
delay between the telescopes. Obviously, any intrinsic star 
variability will be unambiguously identified, as it does not 
produce an intensity geometrical pattern. 

11.3. From space 

The main advantages of a space mission would be the bet- 
ter photometric precision and especially the much better 
spatial resolution, allowing to seriously reduce blending 
problems. As the intensity scintillation due to atmosphere 
is not critical at the level of 1%, the cancellation of this 
specific scintillation in a spatial observatory will be in- 
teresting only for very weak scintillation searches. Among 
the operational spatial observatories and the planned mis- 
sions, the HST observatory and COROT mission have 
been examined: 

— COROT mission (Baglin et al. 2000) plans to monitor 
stars with a short time-scale sampling, but 1% preci- 
sion photometry on My = 20 stars in 10 s exposures 
is definitely unaccessible due to the modest size of the 
telescope. 

— The Space Telescope Imaging Spectrograph instru- 
ment (STIS) of the Hubble Space Telescope (HST) is 
the only instrument that is able to offer a cycle time 
of ~ 20 s for time-resolved spectroscopy or photome- 
try; but this mode is only available on a sub-array of 
the CCD (Profhtt et al. 2002). Regrettably, only one 
(very) wide band filter (550 — 1100 nm) can be used to 
get ~ 5% photometric precision in 10 s exposure time 
on a My = 20.5 BOV star. This very wide band filter 
would already limit the temporal coherence, thus de- 
creasing the diffractive scintillation contrast. The otlic^r 
filters have too narrow passbands to enable short expo- 
sures. Nevertheless, the HST angular resolution would 
be extremely valuable for monitoring M31 fields, and 
there could be an interesting opportunity to further in- 
vestigate. The other HST instruments have prohibitive 
overhead times. 

12. Data flow and analysis 

Data flow and analysis problems can be considered as the 
result of an hybrid of the EROS survey - for the mas- 
sive photometric reduction - and the VIRGO gravitational 
wave experiment. For example, robust flltering developed 
in VIRGO for burst searches (Arnaud et al. 2003) should 



be appropriate to the search for an oscillating signal within 
a definite short period. The time scales of the data flow 
and of the signal are just multiplied by five orders of mag- 
nitude with respect to VIRGO. 

13. Studies to be done, developments and 
perspectives 

A simulation of the Galactic H2 gas distribution, and - 
more important - of the intra-cloud turbulence is clearly 
needed to produce more quantitative predictions on the 
optical depth, on the signal shape and on the temporal 
power spectrum. 

Examination of the inverse problem will lead to more 
precise ideas on the information about the screen and the 
source that could be extracted from a statistical analy- 
sis based on the temporal power spectra. For example, 
as there is a clear connection between the diffractive pat- 
tern, the fringe contrast and the source size, the scintil- 
lation process may also be useful to improve the stellar 
radii knowledge, and to constrain supernovae and quasar 
dimensions. 

It has been amply demonstrated that the angular size 
of the source is the critical parameter for the fringe con- 
trast. Observations with a very large telescope would make 
possible the monitoring of stars much more distant than 
M31, thus providing more contrasted diffractive intensity 
scintillation effects. 

Existing data sets may allow the extraction of inter- 
esting constraints on the refractive scintillation mode. For 
instance EROSl experiment produced 20 000 photomet- 
ric measurements of ~ 100 000 stars in the LMC bar 
(Aubourg et al. 1995). 12 minutes blue and 8 minutes red 
exposiircs were alternatively taken. As refractive scintil- 
lation mode is much less sensitive to the source size than 
diffractive scintillation mode, the sources monitored by 
EROSl were likely to scintillate. As this refractive mode 
could produce long time scale variations (longer than a 
few minutes), the EROSl data set potentially contains 
interesting events. 

As a final remark, one should also consider that the 
systematic survey of thousands of extra-galactic stars at 
high frequency (^ 0.1 Hz) as proposed in this paper may 
give rise to surprises. 

14. Conclusions 

The minimal condition to observe an extra-galactic star 
scintillating through a Galactic molecular cloud is the ex- 
istence of a non vanishing second order derivative of the 
optical path in the transverse plane. To fulfil this con- 
dition in the optical domain, stochastic column density 
fluctuations of 10^® molecules/cm^ on the lO'' — 10^ km 
length scale should take place, that produce optical path 
fluctuations of a fraction of a wavelength. It follows from 
existing models of molecular extended objects that this 
corresponds to column density relative fluctuations of a 
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few ppm per 10'^ — lO"* km of transverse distance, of the 
same order than the average gradient. 

Observations show that structuration of matter is 
present at ah scales, and certainly do not refute the 
eventuality of stochastic fluctuations producing diffractive 
scintillation. 

In this paper, we showed that there is an observational 
opportunity which results from the subtle compromise be- 
tween the arm-lcvcr of interference patterns due to hypo- 
thetic diffusive objects in the Milky- Way and the size of 
the extra-galactic stars. 

The hardware and software techniques required for 
such observations are currently available. Tests are un- 
der way to validate some of the concepts described here. 
If no technical obstacle arises, there is a true opportunity 
to investigate such effects. 
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